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Abstract 

Although  carrier-phase  observations  are  thought  to  drive  the  accuracy  of  time  and 
frequency  transfer  accuracy  obtainable  using  geodetic  GPS  receivers,  there  has  been  evidence 
of  how  colored  noise  on  the  pseudorange  observations  reduces  the  effectiveness  of  the 
technique.  We  present  a  new  parameterization  of  the  GPS  pseudorange  and  carrier-phase 
observations,  called  the  decoupled-clock  model,  where  a  pseudorange  bias  is  estimated  for  each 
clock  in  the  system  as  a  white  noise  parameter.  The  new  parameterization  separates  the 
pseudorange  observation  colored  noise  from  the  carrier-phase  parameters  (ambiguities  and 
clock),  and  also  isolates  the  ambiguities  as  integer-valued  parameters.  While  ambiguous,  the 
carrier-phase  clock  estimates  yield  continuous  frequency  estimates  at  carrier-phase  precision 
over  carrier-phase  continuous  periods.  The  obtained  frequencies  are  consistent  with  those 
implied  by  the  International  GNSS  Service  Rapid  clock  products  without  solution-boundary 
discontinuities.  The  supporting  rationale  for  the  decoupled-clock  model  is  presented  from 
recently  published  literature  and  actual  processing  examples.  The  carrier-phase  clock  solutions 
are  compared  to  the  IGS  Rapid  solutions  and  common  clock  parameterization.  Combination 
methods  are  considered,  including  the  re-introduction  of  observed  pseudoranges  through 
functions  of  the  estimated  pseudorange  biases,  with  the  goal  of  bridging  losses  of  lock  on  the 
constellations,  and  removal  of  the  ambiguity  in  the  carrier-phase  clocks  to  yield  time  transfer 
capability. 
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INTRODUCTION 

The  use  of  geodetic-type  GNSS  reception  hardware  and  processing  algorithms  is  growing  for  time  and 
frequency  transfer  between  atomic  frequency  standards  over  long  distances.  Reference  [1]  provided  a 
very  good  review  of  geodetic  time  and  frequency  transfer.  It  is  now  recognized  that:  a)  the  pseudorange 
measurement  errors,  whether  with  white  or  colored  noise  and  regardless  of  their  source,  are  responsible 
for  the  solution-boundary  discontinuities  [2]  and  b)  the  analysis  of  these  solution-boundary  discontinuities 
gives  an  indication  of  the  quality  of  individual  station  installations  and  may  help  in  identifying  corrective 
actions.  Using  best  practices  for  the  control  of  the  environment  of  the  GNSS  reception  hardware  is  no 
doubt  the  most  effective  method  for  eliminating  most  sources  of  errors  on  pseudoranges,  e.g.  multipath, 
near-field  reflections,  or  temperature  dependencies.  However,  for  various  reasons,  cost  being  the 
principal  one,  these  environment  control  methodologies  cannot  always  be  implemented  in  remote 
locations,  or  for  large  networks.  Also,  from  an  operational  standpoint,  unless  the  reception  equipment 
itself  simultaneously  loses  lock  on  all  satellite  signals,  there  is  a  compelling  argument  not  to  introduce 
discontinuities  in  processing  over  intervals  where  the  hardware  was  able  to  maintain  phase  connection. 

Since  2005,  when  [1]  predicted  the  development  of  tools/software  to  remove  the  burden  of  complex 
geodetic  computations  for  non-geodesists,  many  such  tools  have  appeared  from  various  groups  in  the 
community.  Most  of  these  tools  take  advantage  of  the  International  GNSS  Service  (IGS)  satellite  clock 
and  orbit  products  [3]  in  single-station  processing  and  all  have  two  principal  objectives  in  common:  the 
removal  or  attenuation  of  discontinuities  in  the  clock  series  that  arise  at  solution  time  boundaries  and, 
more  recently,  its  corollary,  the  resolution  of  the  integer  nature  of  undifferenced  carrier-phase 
ambiguities.  To  achieve  this  goal,  the  newly  developed  tools  use  a  range  of  methods:  processing  over 
multiple  days  [4]  and  overlapping  solutions  [5,6],  completely  disregarding  the  pseudoranges  in  carrier- 
phase-only  processing  [7],  resolving  single-differenced  floating  ambiguity  estimates  to  integers  [8],  or 
connecting  the  carrier-phase  ambiguities  between  adjacent  solutions  [9].  Also,  GNSS  frequency  transfer 
results  are  now  combined  with  collocated  independent  time  transfer  techniques  such  as  Two-Way 
Satellite  Time  and  Frequency  Transfer  (TWSTFT)  [10,11]. 

The  issue  of  pseudorange  errors  affecting  clock  estimates  in  GNSS  processing  prompted  us  to  take  a 
closer  look  at  the  fundamental  observation  parameterization  used  in  the  clock  estimation  process.  In  the 
subsequent  sections,  we  will  describe  the  new  decoupled-clock  observation  model  and  the  rationale  that 
led  to  its  development.  We  will  then  test  it  in  a  network-mode  station  and  with  a  satellite  clock 
estimation  least-squares  filter  developed  at  NRCan  for  the  computation  of  GPS  wide-area  corrections. 


GPS  OBSERVATION  MODELING 

The  decoupled-clock  observation  model  has  been  introduced  in  [12,13].  We  will  give  a  summary  of  the 
model  and  the  supporting  rationale. 

Common-Clock  Model 

Dual-frequency  GPS  pseudorange  and  carrier-phase  observations  are  expressed  as  follows: 
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Px  =  p  +  T  +  I  +  c[dtr  -dt')+  b'Pl  - b'Pl  +  sPl 
P ^  =  p  +  T  +  q  I  +  c(dt  —  dt  ^  +  bp~,  —  b p 2  +  £  P2 
A1(Ol+Nl)  =  L1  =  p  +  T  -  I  +  c(dtr -dt')+brn-b‘u+£n 
2  +  N 2)  =  L2  =  p  +  T  -q2I  +  c{dt'  -dt')+b'L2  - b'L1  +£L2 


where  Pi  are  pseudoranges  and  Li  are  carrier  phases  on  each  of  the  GPS  base  frequencies  -  here,  for 
simplicity,  we  do  not  distinguish  between  CA-code  or  P-code  pseudoranges  ([12]  provides  important 
details  on  the  different  pseudorange  types),  p  is  the  geometric  propagation  delay,  T  is  the  propagation 
delay  through  the  troposphere,  /  the  propagation  delay  on  the  L 1  frequency  due  to  the  ionosphere,  q2  the 
conversion  factor  of  ionospheric  propagation  delay  for  the  L2  frequency,  c  the  speed  of  light  in  vacuum, 
dt'  and  dt '  respectively  the  receiver  and  transmitter  clock  offset  with  respect  to  a  common  hypothetical 
reference,  and  the  s  parameters  express  measurement  errors.  The  b*  parameters  are  receiver-  and 
transmitter-specific  bias  parameters  for  each  observable  type.  Dual-frequency  processing  algorithms 
eliminate  the  ionospheric  propagation  delay  to  the  first  order  by  respectively  combining  the  pseudoranges 
and  carrier  phases  on  the  LI  and  L2  transmitted  frequencies  to  generate  the  ionosphere-free  observables: 

_  772P1  -60 2P2 

3_  772  -  602  (Eq.  2) 

—  p  +  T  +  c(dt  —  dt  )  +  b P2  —  b P2  +  s P2 


112Lx  -602L2 

3_  772  -  602  (Eq.  3) 

—  p  +  T  +  c(dt  —  dt  bP2  —  bP2  —  A2N2  +  £ 

where  the  receiver  and  transmitter  clock  parameters  are  the  same  for  pseudorange  and  carrier  phase.  The 
observable  biases  are  usually  considered  negligible  or  as  constants  over  time  absorbed  in  the  clock  offset 
parameters.  Yet  [2]  has  shown  the  pseudorange  biases  are  not  always  negligible  and  more  importantly 
not  always  constant.  When  not  properly  accounted  for  in  the  parameterization,  estimation  processes  will 
distribute  these  time-varying  biases  into  all  estimated  parameters,  primarily  to  those  most  mathematically 
correlated:  clocks  and  ambiguities.  The  direct  effect  is  the  solution-boundary  (usually  daily) 
discontinuities  for  batch  processing,  as  those  observed  in  the  IGS  clock  products  [1],  or  divergence  or 
frequency  errors  in  continuous  filtering  processes.  These  effects  also  reduce  the  ability  of  undifferenced 
data  processing  algorithms  to  resolve  integer  ambiguities.  The  effect  of  the  biases  are  shown  in  Figure  1, 
where  three  relative  clock  solutions  between  two  IGS  stations  (YELL  and  AMC2)  equipped  with 
hydrogen  masers  (HM)  are  displayed:  1)  a  simple  pseudo-range  solution,  2)  a  continuous  least-squares 
filter  using  the  above  common-clock  parameterization  (Eq.  2,3),  and  3)  the  IGS  Rapid  (IGR)  clock 
solution.  The  IGR  solution-boundary  discontinuities  are  present  and  their  sizes  correlate  well  with  the 
amplitude  of  the  long-term  variations  of  the  pseudorange  solution.  Clock  estimates  from  the  filter 
solution  are  continuous,  but  show  frequency  errors  when  compared  to  the  IGR  (positive  vs.  negative 
slopes  and  vice  versa).  Note  that,  since  within  one  daily  batch  the  IGR  solutions  are  combinations  of 
ambiguity  fixed  solutions,  these  should  represent  the  clock  as  ‘seen’  by  the  carrier  phases.  We  will 
hereafter  label  this  parameterization  and  any  results  thereof  “common-clock.” 
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YELL-AMC2  clock  error 
(Jan07-Janl3f2007;  common  Unear  lit  removed) 


Dqr  of  Week  (GES 1409) 


Figure  1.  Three  clock  solutions  for  the  YELL-AMC2  station  pair  involving  two  HM 
(from  [12]). 


Decoupled-Clock  Model 


Taking  into  account  the  biases  in  the  parameterization  leads  to  a  new  system  of  equations: 


P3  =  p(x  , x  )  +  zpd  ■  mapT  +  c\dt'a  -  dt'L3 )+  B'n  -  B'p 


L3  =  p{x  ,x  )  +  zpd  ■  mapT  +  c\dt'n  -  dt'Li 


+  vf 


A3N3  +vu 


(Eq.  4) 


where  overbars  now  express  estimated  parameters,  p()  is  the  distance  between  transmitter  and  receiver 

-*r  -*t 

instantaneous  positions  (x  ,x  )  accounting  for  earth  displacements  and  including  relativistic  corrections, 
mapT  is  the  function  mapping  the  local  receiver  zenith  tropopheric  delay  zpd  to  the  slant  delay  T,  dt' L3  and 
dt‘u  are  the  receiver  and  transmitter  carrier-phase  clock  offsets,  B‘ Pi  and  B' P3  are  the  receiver  and 
transmitter  pseudorange  biases  (or  decoupled  “pseudorange  clock  offsets”),  13N3  is  the  L  ;  carrier-phase 
ambiguity  multiplied  by  its  wavelength,  and  vP3  and  vL3  the  pseudorange  and  carrier-phase  observation 
residual  errors.  Because  X3  is  ~6mm,  too  small  to  allow  efficient  integer  ambiguity  resolution  if  at  all 
possible,  we  add  the  Melbourne-Wubenna  observation  combination  [14,15]: 


A4  =  L4-P6  =(77^  -60Z2)/17  ~(77Pl  -6(TP,)/137 
=  BA4  -  Ba4  —  A4  N4  +  vA4 


(Eq.  5) 
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composed  of  the  wide-lane  phase  combination  L4  and  the  pseudorange  narrow-lane  combination  P6.  The 
)-4  wavelength  of  ~86  cm  allows  easier  integer  resolution  for  N4  which  in  turn,  through  parameter 
substitution,  gives  a  wavelength  of  -10.7  cm  for  N3,  better  suited  for  integer  resolution  techniques. 

As  it  stands,  the  system  of  equations  formed  by  Eqs.  4  and  5  is  singular,  the  biases  and  clock  terms  being 
linearly  dependent.  By  decoupling  the  pseudorange  and  carrier-phase  clock  terms,  we  effectively 
removed  the  datum  provided  by  the  pseudorange  measurements  which  allowed  ambiguity  estimation  in  a 
common-clock  parameterization.  We  resolve  this  rank  deficiency  by  providing  an  arbitrary  ambiguity 
datum  and  this  will  yield  ambiguous  carrier-phase  clock  estimates.  To  provide  a  set  of  minimum- 
constraint,  solvable  normal  equations,  one  ambiguity  per  carrier-phase  clock  less  one,  must  be  fixed  to  an 
arbitrary  value.  Additionally,  one  set  of  carrier-phase  clock,  pseudorange  bias,  and  wide-lane  bias 
parameters  must  be  fixed  as  the  network  datum.  This  is  similar  to  network  processing  using  the  common- 
clock  parameterization  where  one  clock  must  be  fixed.  To  maintain  the  continuity  of  all  parameters,  this 
minimum  set  of  fixed  ambiguities  and  parameters  must  be  carried  over  time,  as  satellites  rise  and  set  over 
receiving  stations.  This  requires  ambiguity  bookkeeping  algorithms  for  on-the-fly  selection  of  integer 
resolved  ambiguities  to  replace  datum  ambiguities  disappearing  from  the  system. 

The  decoupled-clock  parameterization  can  be  exploited  in  both  network  mode,  where  all  transmitter  and 
receiver  clocks  are  solved  for  (less  one  datum  set  as  we  have  already  seen),  as  well  as  in  a  single  receiver 
Precise  Point  Positioning  (PPP)  mode  [16,17],  as  was  shown  in  [13].  The  requirement  here  is  a 
decoupled-clock  network  solution  for  the  transmitter  constellation,  e.g.,  carrier-phase  clocks,  pseudorange 
and  wide-lane  biases,  and  fixing  a  single  ambiguity  parameter  to  provide  the  ambiguity  datum.  It  was 
also  shown  in  [131  that  ambiguity  resolution  techniques  could  effectively  be  applied  to  the  decoupled- 
clock  parameterization. 

The  decoupled-clock  model  was  implemented  in  a  sequential  least-squares  filter  which  was  originally 
developed  for  the  generation  of  GPS  wide-area  corrections  [18].  The  filter  accepts  observations  either 
from  real-time  data  streams  or  from  flat  files  in  the  estimation  of  satellite  and  station  clocks,  local  station 
troposphere  delay,  and  carrier-phase  ambiguities,  to  which  we  add  the  pseudorange  and  wide-lane  biases 
when  using  the  decoupled-clock  model.  Stations  and  satellites  positions  are  held  fixed  and  a 
comprehensive  set  of  conventions  and  physical  models  are  adopted  and  used.  Since  the  least-squares 
filter  is  continuous,  it  does  not  have  solution-boundary  issues.  Using  the  common-clock  model,  the  effect 
of  pseudorange  measurement  errors  is  that  the  ambiguities  become  so  much  constrained  to  slightly 
erroneous  values  through  accumulation  of  weight  over  time,  that  the  system  starts  diverging  after  about  a 
month  or  more.  We  initially  resolved  the  divergence  issue  by  adding  small  process  noise  to  the  ambiguity 
parameters.  This,  however,  is  not  entirely  satisfactory,  as  the  pseudorange  errors  contaminate  the  clock 
solutions  through  the  floating  ambiguities,  as  was  shown  in  Figure  1  and  will  be  seen  later. 


EXPERIMENTAL  ASSESSMENT 

For  the  assessment  of  the  new  observation  model,  the  least-squares  filter  was  used  to  process  the 
following  dataset.  Station  coordinates  were  held  fixed  to  their  IGS05  coordinates  of  epoch  and  satellites 
to  the  IGR  orbits.  The  GPS  PI  and  P2  pseudorange  observation  pair  was  used  in  the  ionosphere-free 
combination,  consistent  with  the  IGS  clock  product  convention.  Receivers  producing  only  CA 
pseudoranges  are  corrected  for  the  Cl-Pl  biases  [19].  The  IGS  convention  for  satellite  and  station 
absolute  antenna  phase  center  offsets  and  azimuth-elevation  variations  definition  was  used  as  well  [20]. 

Two  solutions  were  generated  with  our  processing  engine:  1)  using  the  common-clock  model  with 
introduction  of  process  noise  in  the  ambiguities  and  2)  using  the  decoupled-clock  model. 
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Data  Set 

We  selected  a  data  set  from  the  IGS  archive  of  daily  30-second  RINEX  format  [21]  files  to  form  a 
reception  network  of  29  stations  distributed  around  the  globe,  spanning  35  days  in  early  2008,  specifically 
from  MJD  54499  to  54533  inclusive,  to  which  we  added  1 1  stations  of  interest  for  time  and  frequency 
transfer  purposes.  The  selection  criteria  for  these  11  stations  were:  1)  timing  laboratories  operating 
geodetic  GNSS  receivers,  2)  stations  using  atomic  frequency  standards,  and  3)  stations  presenting  a  good 
reception  record  for  the  target  period.  The  selection  criteria  for  the  29  other  stations  were  simply  to 
ensure  continuous  tracking  of  the  satellite  constellation.  This  is  of  utmost  importance  for  maintaining  the 
ambiguity  datum  over  time.  Table  1  provides  the  list  of  selected  stations  and  their  GNSS  reception 
equipment  as  reported  by  the  logs  maintained  at  IGS.  IGR  clock  products,  aligned  to  the  IGR  timescale, 
were  also  collected  for  comparison. 


Pseudorange  Measurement  Errors 

An  example  of  clock  results  contaminated  by  pseudorange  measurement  errors  is  shown  in  Figure  2, 
where  the  ambiguous  phase  clock  solution  is  compared  to  the  common-clock  model  and  IGR  clocks  for 
the  YELL  and  PTBB  station  pair.  Figure  3  shows  the  pseudorange  and  wide-lane  biases  estimated  in  the 
decoupled-clock  model  filter  process  for  the  same  period  and  station  pair.  The  pseudorange  biases  show 
large  variations  between  MJD  54503  and  54510  and  between  MJD  54526  and  54531.  The  response  of 
the  continuous  common-clock  model  processing  is  to  follow-on  average  the  pseudorange  biases  and  that 
of  the  IGR  product  is  larger  day-boundary  discontinuities.  The  carrier-phase  clocks  from  the 
decoupled-clock  model  are  continuous  over  carrier-phase-continuous  periods  and  their  frequency  is 
similar  to  that  of  IGS  within  each  day.  Daily  averages  of  the  estimated  pseudorange  biases  were 
computed  as  well  as  an  average  over  carrier-phase-continuous  periods  (Figure  3).  There  is  a  strong 
correlation  between  the  IGR  product  day-boundary  discontinuities  and  those  implied  by  the  daily 
averages  of  pseudorange  bias.  The  correlation  is  even  more  striking  in  Figure  4,  where  the  carrier-phase 
clocks  were  differenced  with  IGR  clocks  and  then  with  IGR  clocks  corrected  for  the  jumps  implied  by  the 
daily  pseudorange  bias  averages.  Correcting  the  IGR  clocks  day-boundary  discontinuities  using  daily 
averages  of  the  pseudorange  biases  as  estimated  by  the  decoupled-clock  model  brought  the  consistency  of 
both  solutions  at  the  sub-nanosecond  level.  The  carrier-phase  clock  from  the  decoupled-clock  model  is 
largely  immune  to  the  pseudorange  bias  variations,  although  there  are  still  half-daily  oscillations  in  the 
differences  that  can  be  attributed  to  the  carrier-phase  solution,  as  will  be  shown  in  the  stability  analysis. 

Table  2  lists  RMS  differences  of  nine  station  pairs  for  the  decoupled-model  carrier-phase  clocks  and  IGR 
clocks,  both  corrected  and  uncorrected  for  day-breaks  as  explained  above.  Overall,  the  consistency  of  the 
decoupled  carrier-phase  clocks,  the  pseudorange  biases  and  IGR  clocks  is  at  the  200-ps  RMS  level,  or 
better.  One  should  notice  that  the  YELL-PTBB  results  are  an  extreme  case  caused  by  local  conditions 
prevalent  at  the  YELL  station  and  well  documented  in  [1]. 

For  the  other  stations  pairs,  where  the  pseudorange  errors  do  not  have  such  significant  signatures,  the 
improvement  in  RMS  differences  by  correcting  IGR  clocks  for  day-breaks  is  not  as  pronounced.  The 
overall  improvement  is,  on  average,  about  50%,  even  when  excluding  the  two  extreme  cases  (NRC1, 
YELL).  This  is  a  reflection  of  the  internal  consistency  of  the  IGS  day-boundary  discontinuities  and 
decoupled-clock  model  pseudorange  bias  estimates  and  gives  us  confidence  the  new  parameterization  is 
sound. 
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Table  1.  IGS  stations  selected  and  equipment  list  from  MJD  54499  to  54533  as  reported 
at  http://igs.org. 


Station 

Receiver 

Antenna 

Frequency 

Standard 

Notes 

Timing  laboratories 

amc2 

ashtech  Z-XII3T 

aoad/m  t 

h-maser 

BRUS 

ashtech  Z-XII3T 

ash701945b  m 

h-maser  ch1-75 

IENG 

ASHTECH  Z-XII3T 

ash701945c  m 

h-maser 

NRCl 

AOA  SNR- 12  ACT 

aoad/m  t 

h-maser 

NRLl 

ASHTECH  Z-XII3T 

ash701945c  m 

h-maser 

PTBB 

ASHTECH  Z-XII3T 

ash700936e 

cesium 

SNOW  RADOME 

SPTO 

JPS  LEGACY 

aoad/m_t 

h-maser 

OSOD  RADOME 

SYDN 

JPS  E  GGD 

ash701945c_m 

cesium 

TWTF 

ASHTECH  Z-XII3T 

ash701945c_m 

steered  h-maser 

SCIS  RADOME 

usn3 

ASHTECH  Z-XII3T 

aoad/m  t 

h-maser 

wab2 

ASHTECH  Z-XII3T 

ash700936f  c 

h-maser 

SNOW  RADOME 

Other  atomic  FS 

coco 

ASHTECH  UZ-12 

aoad/m  t 

RUBIDIUM 

DARW 

ASHTECH  UZ-12 

LEICA  GRX1200GGPRO 

ash700936d_m 

RUBIDIUM 

DAV 1 

ASHTECH  UZ-12 

ash701945g  m 

RUBIDIUM 

AUST  RADOME 

DRAO 

AOA  BENCHMARK  ACT 

aoad/m  t 

RUBIDIUM 

KOKB 

ASHTECH  UZ-12 

ash701945g  m 

H-MASER 

KOUR 

JPS  LEGACY 

ASH701946.3 

CESIUM 

LPGS 

AOA  BENCHMARK  ACT 

aoad/m  t 

Rubidium 

MASl 

ASHTECH  Z-XII3 

aoad/m  t 

CESIUM 

MIZU 

SEPT  POLARX2 

TPSCR3  GGD 

Cesium 

OUS2 

SEPT  POLARX2 

TPSCR3  GGD 

RUBIDIUM 

SANT 

ASHTECH  UZ-12 

aoad/m  t 

CESIUM 

JPLA  RADOME 

SUTM 

AOA  BENCHMARK  ACT 

aoad/m  t 

RUBIDIUM 

TIDB 

ASHTECH  Z-XII3 

aoad/m  t 

H-MASER 

JPLA  RADOME 

tow2 

ASHTECH  UZ-12 

aoad/m  t 

RUBIDIUM 

AUST  RADOME 

YELL 

AOA  SNR- 12  ACT 

aoad/m  t 

H-MASER 

Other  stations 

CHPI 

ASHTECH  UZ-12 

ash701945c_m 

quartz/internal 

FAAl 

JPS  LEGACY 

ash701945b_m 

quartz/internal 

HYDE 

LEICA  GRX1200GGPRO 

leiat540gg 

UNSPECIFIED 

ISPA 

ASHTECH  UZ-12 

ash701945e  m 

QUARTZ/INTERNAL 

SCIT  RADOME 

KARR 

ASHTECH  UZ-12 

aoad/m  t 

QUARTZ/INTERNAL 

AUST  RADOME 

MALI 

JPS  LEGACY 

ASH701946.3 

UNSPECIFIED 

MKEA 

ASHTECH  Z-XII3 

aoad/m  t 

QUARTZ/INTERNAL 

OHI2 

JPS  E  GGD 

aoad/m  t 

QUARTZ/INTERNAL 

UNSPEC.  RADOME 

ohi3 

TPS  E  GGD 

ASH701941.B 

QUARTZ/INTERNAL 

SNOW  RADOME 

RCMN 

LEICA  GRX1200GGPRO 

leiat540gg 

UNSPECIFIED 

LEIS  RADOME 

SEY 1 

ASHTECH  UZ-12 

ash701945b  m 

QUARTZ/INTERNAL 

UNSPEC.  RADOME 

TASH 

Unspecified 

ULAB 

AOA  snr-8000  act 

aoad/m  t 

UNSPECIFIED 
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54499  54506  54513  54520  54527  54534 

Figure  2.  Thirty-five-day  decoupled-clock  solution  for  the  YELL-PTBB  station  pair. 
Common-clock  model  and  IGR  solutions  included  for  comparison.  All  solutions  have  a 
common  frequency  offset  and  drift  removed  and  are  arbitrarily  offset  in  phase  for  display. 
Daily  average  of  pseudorange  biases  included. 


Table  2.  RMS  of  differences  between  decoupled-clock  model  carrier-phase  clocks  and 
IGR  corrected  and  uncorrected  for  day-breaks.  The  “Partial”  row  excludes  pairs  with 


NRC1  and  YELL. 


Station  pair 

#  diff 

PhaseClock-IGS 
RMS  (ns) 

PhaseClock  -  IGScorr. 
RMS  (ns) 

Improv. 

% 

AMC2-PTBB 

7591 

0.142 

0.111 

22.0 

BRUS-PTBB 

9375 

0.195 

0.064 

67.1 

IENG-PTBB 

9361 

0.103 

0.066 

35.8 

KOKB-PTBB 

9377 

0.237 

0.135 

43.1 

NRC1-PTBB 

9314 

0.729 

0.197 

73.0 

SPTO-PTBB 

9294 

0.209 

0.111 

46.7 

TWTF-PTBB 

9377 

0.220 

0.129 

41.4 

USN3-PTBB 

8428 

0.263 

0.119 

54.9 

YELL-PTBB 

8789 

1.365 

0.200 

85.4 

Partial 

72167 

0.247 

0.115 

53.4 

Overall 

90270 

0.534 

0.136 

74.5 
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54506 


54513 


54520 


-40 
54534 


Figure  3.  Pseudorange  and  wide-lane  biases  for  the  same  period  and  station  pair.  Daily 
averages  (black)  of  the  estimated  pseudorange  biases  as  well  as  straight  average  (red) 
over  the  three  phase-continuous  periods  are  plotted. 


The  noise  characteristics  of  the  estimated  pseudorange  biases  were  further  analyzed  using  Overlapping 
Allan  deviations,  as  shown  in  Figure  5.  The  Overlapping  Allan  deviations  show  little  differences  in  the 
short  term  and  only  small  differences  occur  over  the  longer  averaging  intervals.  Further  analysis  will  be 
required  to  investigate  their  significance. 
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Figure  4.  Differences  between  carrier-phase  clock  and  IGR  (blue)  and  IGR  corrected  for 
the  day-breaks  implied  by  the  daily  averages  of  pseudorange  biases  (red). 
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Figure  5.  Overlapping  Allan  deviation  of  pseudorange  bias  estimates  for  all  station  pairs, 
datum  jumps  removed.  Tau'1  noise  line  shown. 


Stability  Analysis 

We  performed  a  stability  analysis  of  the  carrier-phase  clock,  computing  their  Overlapping  Allan 
deviations.  For  each  of  the  baselines  listed  in  Table  3,  we  selected  the  stated  periods  within  the  35  days 
where  there  were  no  clock  resets,  no  ambiguity  datum  resets,  and  no  missing  days  in  the  IGS  clock 
solutions.  Frequency  series  were  generated  from  the  following  clock  solutions:  1)  the  decoupled-clock 
model  carrier-phase  clock,  2)  the  carrier-phase  clocks  where  we  impose  the  day-breaks  implied  by  the 
daily  averages  of  estimated  pseudorange  biases,  3)  the  common  clock  solution,  4)  the  IGR  clock  solution, 
and  5)  the  IGR  clock  model,  excluding  those  intervals  crossing  day-boundaries.  The  last  case  is  not 
realistic  and  here  considered  an  ideal  case,  i.e.  IGS  quality  clocks  without  day-boundary  discontinuities 
and,  thus,  uncontaminated  by  pseudorange  errors.  Data  gaps  with  respect  to  the  nominal  solution 
intervals  (300  seconds  for  the  IGR  and  30  seconds  for  the  other  series)  were  considered  as  missing  data. 
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Table  3.  Selected  station  pairs  and  associated  periods  for  stability  analysis. 


Station  pair 

Period 

#days 

AMC2-BRUS 

54517.69-54529.53 

11.84 

IENG-BRUS 

54519.78  -  54529.53 

9.74 

KOKB-BRUS 

54501.33  -54515.77 

14.44 

PTBB-BRUS 

54501.33  -54529.53 

12.66 

SPTO-BRUS 

54515.78  -  54529.53 

13.75 

TWTF-BRUS 

54515.78  -  54529.53 

13.75 

USN3-BRUS 

54515.78  -  54529.53 

13.75 

YELL-BRUS 

54509.00-54515.77 

6.77 

Figure  6  shows  the  Overlapping  Allan  deviations  for  the  YELL-BRUS  and  AMC2-BRUS  station  pairs. 
For  the  YELL-BRUS  pair,  affected  most  by  YELL  pseudorange  biases,  the  improvement  obtained  by 
removing  the  day-boundary  discontinuities  from  the  IGR  clocks  is  almost  an  order  of  magnitude.  The 
carrier-phase  clock  stability,  while  not  as  good  as  the  ideal  case,  is  better  than  that  obtained  from  the 
common  clock  model.  Introducing  day-breaks  in  the  otherwise  continuous  carrier-phase  clocks  brings  its 
stability  to  the  level  of  the  IGS  solution.  For  the  AMC2-BRUS  pair,  the  common  clock  and  decoupled- 
clock  model  solutions  have  similar  stability,  a  little  worse  than  the  IGR.  It  appears  the  day-boundary 
discontinuities  have  little  effect  on  that  baseline.  Worthy  of  notice  is  the  small  rise  in  the  Allan  deviations 
for  the  decoupled-clock  model  around  a  6-hour  interval.  This  confirms  that  the  small  half-day 
oscillations  seen  in  Figure  4  come  from  the  decoupled-clock  model.  This  will  be  the  subject  of  further 
analysis.  Tables  4  and  5  respectively  show  the  stability  at  1-day  and  1-hour  for  all  the  station  pairs 
studied.  For  all  station  pairs,  the  carrier-phase  clocks  have  improved  long-term  stability  over  the  IGR  and 
are  close  to  the  ideal  case  represented  by  IGR  without  day-boundaries.  For  most  cases,  it  is  as  good  or 
better  than  the  common-clock  results. 


Table  4.  Overlapping  Allan  deviation  at  1  day  for  the  frequency  series  and  station  pairs 
(units  10'15). 


Station  pair 

IGS  without 
day-breaks 

Carrier-phase 

clock 

IGS  with 
day-breaks 

Carrier-phase 
with  breaks 

Common 

clock 

AMC2-BRUS 

2.5 

2.5 

2.8 

2.8 

2.5 

IENG-BRUS 

2.7 

2.6 

2.7 

2.7 

2.5 

KOKB-BRUS 

2.4 

2.4 

2.7 

2.8 

2.4 

PTBB-BRUS 

12.1 

12.2 

13.0 

12.6 

12.3 

SPTO-BRUS 

4.0 

4.1 

5.8 

5.9 

4.3 

TWTF-BRUS 

3.7 

3.6 

3.8 

3.5 

3.9 

USN3-BRUS 

2.3 

2.3 

2.2 

2.2 

1.9 

YELL-BRUS 

3.9 

4.1 

21.4 

22.2 

6.9 
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Table  5.  Overlapping  Allan  deviation  at  1  hour  for  the  frequency  series  and  station  pairs 
(units  10'14). 


Station  pair 

IGS  without 
day-breaks 

Carrier-phase 

clock 

IGS  with 
day-breaks 

Carrier-phase 
with  breaks 

Common 

clock 

AMC2-BRUS 

1.2 

1.9 

1.3 

1.9 

1.9 

IENG-BRUS 

1.4 

1.4 

1.5 

1.5 

1.9 

KOKB-BRUS 

1.3 

2.5 

1.6 

2.6 

2.1 

PTBB-BRUS 

5.6 

5.8 

5.8 

5.8 

5.8 

SPTO-BRUS 

1.3 

1.2 

1.8 

1.8 

1.4 

TWTF-BRUS 

1.2 

2.0 

1.5 

2.0 

2.0 

USN3-BRUS 

1.2 

1.8 

1.5 

1.9 

1.7 

YELL-BRUS 

1.0 

1.6 

8.7 

9.0 

3.4 

CONCLUSION 

The  decoupled-clock  model,  a  new  parameterization  of  the  GPS  carrier-phase  and  pseudorange 
observations,  and  the  rationale  behind  its  development  were  described.  The  objective  of  the  new 
parameterization  is  to  remove  the  influence  of  pseudorange  errors  on  the  receiver  clock  estimates,  which 
in  batch  solutions  take  the  form  of  batch-boundary  discontinuities  which  are  too  large  to  be  accounted  by 
noise  characteristics  of  the  atomic  frequency  standard  to  which  the  receiver  is  connected.  The  decoupled- 
clock  model  was  implemented  in  a  least-squares  filter  designed  for  network-mode  station  and  satellite 
clock  determinations  for  the  generation  of  wide-area  GPS  corrections.  Ambiguous  carrier-phase  clocks 
and  pseudorange  biases  were  obtained  by  processing  35  days  of  tracking  data  from  a  network  of  39 
globally  distributed  IGS  stations,  1 1  of  which  are  located  at  timing  laboratories  and  several  others  are 
connected  to  atomic  frequency  standards.  The  same  data  set  was  processed  using  the  common 
parameterization  where  pseudorange  biases  are  not  estimated  and  with  process  noise  applied  to  the 
carrier-phase  ambiguities  estimation  to  prevent  divergence. 

Daily  averages  of  the  estimated  pseudorange  biases  are  consistent  with  the  day-boundary  discontinuities 
present  in  the  IGR  clock  solutions.  Using  these  averages  to  correct  the  IGR  clock  solutions  brought  their 
consistency  with  the  continuous  decoupled-clock  model  carrier-phase  clock  solutions  at  the  sub-ns  level 
even  for  the  most  affected  stations  YELL  and  NRC1.  Lor  the  other  stations  where  pseudorange 
contamination  is  not  readily  apparent,  the  clock  solution  consistency  has  improved  by  50%  over  all  the 
stations  analyzed,  giving  us  confidence  in  the  new  clock  parameterization. 

The  stability  of  the  carrier-phase  clock  solutions  were  analyzed  in  terms  of  Overlapping  Allan  deviations 
for  the  IGR  clocks  both  with  and  without  the  intervals  crossing  the  problematic  day-boundaries,  the  later 
being  an  ideal  case.  Lor  stations  most  contaminated  by  pseudorange  biases,  the  improvement  brought  by 
the  decoupled-clock  model  is  significant,  close  to  an  order  of  magnitude  in  case  of  station  YELL.  Lor 
other  stations,  the  improvement  in  stability  is  less  significant. 
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Figure  6.  Overlapping  Allan  deviation  for  frequency  series  from  YELL-BRUS  and 
AMC2-BRUS  station  pairs.  Tau'12  noise  line  shown. 
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The  objective  of  separating  pseudorange  errors  within  clock  and  ambiguity  estimation  is  largely  met  by 
the  decoupled-clock  model,  although  for  the  most  affected  stations,  there  are  still  half-daily  oscillations  in 
the  carrier-phase  clocks  that  need  further  investigations.  At  its  current  stage  of  development,  since  the 
obtained  clocks  are  ambiguous,  the  new  clock  parameterization  yields  frequency  transfer  accuracy  that  is 
consistent  with  the  one  implied  by  the  IGR  clocks  without  day-boundary  discontinuities,  for  as  long  as 
both  receivers  maintain  carrier-phase  lock  on  the  satellite  constellation.  The  decoupled-clock  ambiguity 
can  be  removed  by  combining  with  collocated  independent  time  transfer  methods,  e.g.  TWSTFT. 

There  is  a  side  benefit  that  far  exceeds  just  meeting  the  objective  we  have  set.  The  new  model  explicitly 
reveals,  on  average,  the  noise  structure  of  the  pseudoranges,  rather  than  seeing  their  effects  compounded 
into  the  clock  results  as  solution-boundary  discontinuities.  This  opens  areas  of  investigation  and  future 
development.  Analysis  of  the  noise  properties  on  the  pseudorange  biases  would  allow  the  selection  of  the 
most  effective  function  thereof  to  reintroduce  the  pseudorange  datum  in  the  results.  For  example,  if  the 
pseudorange  biases  for  a  particular  station  pair  are  white  phase  noise,  simply  adding  their  average  over 
carrier-phase  continuous  periods  to  the  carrier-phase  clocks  would  allow  connecting  periods  between 
carrier-phase  losses  of  lock,  thus  extending  the  frequency  transfer  capability  over  very  long  continuous 
periods.  Furthermore,  if  the  ambiguity  datum  can  be  resolved  to  the  correct  integer,  the  carrier-phase 
ambiguity  could  be  removed,  yielding  continuous  carrier-phase  quality  GPS  time  transfer. 
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